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Abstract 



We review the method of stochastic quantization for a scalar field theory. We first give a 
brief survey for the case of self-interacting scalar fields, implementing the stochastic perturbation 
theory up to the one-loop level. The divergences therein are taken care of by employing the usual 
prescription of the stochastic regularization, introducing a colored random noise in the Einstein 
relations. We then extend this formalism to the case where we assume a Langevin equation with 
a memory kernel. We have shown that, if we also maintain the Einstein's relations with a colored 
noise, there is convergence to a non-regularized theory. 
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1 Introduction 



In the last century Parisi and Wu introduced the stochastic quantization [1]. This new quan- 
tization method differs from the others, the canonical and the path integral field quantization, 
based in the Hamiltonian and the Lagrangian, respectively, in many aspects. The method starts 
from a classical equation of motion, but not from Hamiltonian or Lagrangian, and consequently 
can be used to quantize dynamical systems without canonical formalism. Furthermore, it is useful 
in situations where the others methods lead to difficult problems and can bring us new important 
results. As stressed by Rumpf [2], since the stochastic quantization is quite different from the 
other quantization methods, it can reveal new structural elements of the classical theory which 
so far have gone unnoticed. The main idea of the stochastic quantization is that a c?-dimensional 
quantum system is equivalent to a (d + l)-dimensional classical system with random fluctuations. 
Some of the most important papers in the subject can be found in Ref. [3]. A brief introduction 
to the stochastic quantization can be found in the Ref. [4] and Ref. [5]. See also the Ref. [6]. 

In this paper we would like to discuss the stochastic quantization method for the (A(/?^)d scalar 
theory using a Langevin equation with a memory kernel and a colored noise. We inquiry which 
kind of theory appears in the asymptotic limit of this non-Markovian process at one-loop level. 
Does this stationary, Gaussian, non-Markovian theory reaches the equilibrium with the structure 
of its ultraviolet divergences under control? 

The method of stochastic quantization can be summarized by the following steps. First, start- 
ing from a field defined in Minkowski spacetime, after analytic continuation to imaginary time, 
the Euclidean counterpart, i.e., the field fiving in an Euclidean space, is obtained. Second, it is 
introduced a monotonically crescent Markov parameter, called in the literature "fictitious time" 
and also a random noise field r]{T,x), which simulates the coupling between the classical system 
and a heat reservoir. It is assumed that the fields defined at the beginning in a ci- dimensional Eu- 
clidean space also depends on the Markov parameter, therefore the field and a random noise field 
are defined in a {d+ l)-dimensional manifold. One starts with the system out of equilibrium at an 
arbitrary initial state. It is then forced into equilibrium assuming that its evolution is governed 
by a Markovian Langevin equation with a white random noise field. In fact, this evolution is de- 
scribed by a process which is stationary, Gaussian and Markovian. Finally, the n-point correlation 
functions of the theory in the {d + l)-dimensional space are defined by performing averages over 
the random noise field with a Gaussian distribution, that is, performing the stochastic averages 
{ (P{ti, Xi)(p{T2, X2) . ■ .(p{Tn, Xn) The n-poiut Schwinger functions of the Euclidean d-dimensional 
theory are obtained evaluating these rz-point stochastic averages {(fi{Ti,Xi)ip{T2,X2)...(p{Tn,Xn))n 
when the Markov parameter goes to infinity (r — > 00), and the equilibrium is reached. This 
can be proved in different ways for the particular case of Euclidean scalar field theory. One can 
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use, for instance, the Fokker-Planck equation [7] [8] associated with the equations describing the 
stochastic dynamic of the system. A diagrammatical technique [9] has also been used to prove 
such equivalence. 

The original method proposed by Parisi and Wu was extended to include process with fermions 
[10] [11] [12]. The first question that appears in this context is if make sense the Brownian problem 
with anticommutating numbers. It can be shown that, for massless fields, there will not be a 
convergence factor after integrating the Markovian Langevin equation. Therefore the equilibrium 
is not reached. One way of avoiding this problem is to introduce a kernel in the Langevin equation 
describing the evolution of two Grassmannian fields. 

Usually, the Parisi- Wu schemes of quantization applied to bosonic and fermionic fields converge 
towards non-regularized theories. In order to obtain finite results, the original stochastic process is 
modified to a non-Markovian process, taking into account a colored noise in the Einstein relations 
[13] [14]. In few words, the idea of the stochastic regularization is to start from an interacting 
theory, then construct Langevin tree-graphs, where each leg of which ends in a regularized noise 
factor. Since it is possible to obtain the loops of the theory by contracting the noise factors, one 
ends up with a theory where every closed loop contains at least some power of the regulator. With 
this modification, one can show that the system converges towards a regularized theory. The next 
step to construct a finite theory is to use for example a minimal-subtraction scheme in which the 
ultraviolet divergent contributions are eliminated. 

However, different types of non-Markovian processes are allowed in principle; for instance, we 
can modify the Langevin equation introducing a memory kernel. To make sure that the generalized 
Langevin equation can also be used as a quantization tool, one must check that the process 
converges in the asymptotic limit and also that converges to the correct equilibrium distribution. 
In this paper we address the following question: which kind of interacting field theory appears 
in the asymptotic limit of this non-Markovian process at one-loop level? We have shown that 
although a system with a stationary, Gaussian, non-Markovian Langevin equation with a memory 
kernel and a colored noise converges to equilibrium, we obtain a non-regularized theory. 

The organization of the paper is the following: in section 2 we briefly discuss the Markovian 
and the non-Markovian Langevin equation in stochastic processes. We discuss in section 3 the 
stochastic quantization for the free scalar theory. The stochastic quantization for the (A(/?^)d self- 
interaction scalar theory with the usual stochastic regularization method is presented in section 
4. In section 5 assuming a Langevin equation with a memory kernel, we study the free and the 
interacting field theory that appears in the asymptotic limit of the non-Markovian process. In 
section 6, we developed a Fokker-Planck approach in order to understand the results obtained in 
the previous section. Conclusions are given in section 7. To simplify the calculations we assume 
the units to be such that h = c = ks = I- 
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2 Markovian and non-Markovian Langevin equation 



Let us briefly discuss the role of the Langevin equation in non- equilibrium statistical mechanics 
[15]. One way to treat the dynamics of non-equilibrium systems is to use the theory of Brownian 
motion. The fundamental equation is the Langevin equation, and it contains both frictional and 
random forces. The equation of motion for a Brownian particle is given by 

mj^^r{t)^-^vit)+5Fit), (1) 

where —^\{t) is the frictional force and SF{t) is the stochastic or fluctuating force, also called 
the random force. In Eq.(l), m is the mass of the Brownian particle immersed in the fluid. Note 
that the total force that acts on each Brownian particle has been partioned into a friction part 
and also a fluctuating part. The effects of the fluctuation force can be summarized by given the 
Einstein relations, namely, 

(5F(t)) = (2) 

and also 

{SF\t) 5F^{t') ) = 2B5ijS{t - t'), (3) 

where B is the measure of the strength of the fluctuating force. The Eq.(2) tell us that the random 
force is zero on average and the delta function in Eq.(3) indicates that there is no correlation 
between impacts in any two distinct time intervals. Since the frictional force depends just on the 
velocity of the particle and not on its earlier values, we are describing a Markovian process. The 
Langevin equation can be solved to give 



v(t) = exp ( - — t ) v(0) + — dsex.p[-—(t- s)] 5F(s). 
\ m J m Jo \ m J 



(4) 



Using the fact that the fluctuating force has a Gaussian distribution defined by its first and seconds 
moments, it is possible to show that the mean squared velocity of each Brownian particle is 



{v'{t) ) = ew[--t] v\0) + — 1 - exp U , (5) 

\ m J qiTi \ \ m J J 

where the brackets (...) in the Eq. (5) denote stochastic averaging. 

There are two points that we would like to stress. First, it is assumed that the fluctuating force 
SF{t) has a Gaussian distribution. Therefore it is possible to relate the strength B of the random 
noise to the magnitude ^ of the friction [B = ^T), where T is the temperature of the bath. This 
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is the well known fluctuation-dissipation theorem. Note that the Langevin equation considered in 
Eq.(l) is called Markovian, since the frictional force {—^"v{t)) at t is proportional to the velocity 
at t and the noise is delta-function correlated or white. Being more specific, the Fourier transform 
of the correlation function of the noise, that is, the spectral density, is independent of frequency. 

The Langevin equation may be generalized to non-Markovian types assuming that the friction 
at time t can depend on the history of the velocity v(s) for times s that are earlier than t [16] 
[17] [18]. In this situation, the friction coefficient must be replaced by a memory function. The 
Langevin equation with a dissipative memory kernel, that we will call a generalized Langevin 
equation, reads 

d /'°° 

m—Y(t)^- dsM(s)v(t- s) + SF(t). (6) 
at Jo 

In the paper we will assume a non-Markovian Langevin equation with a memory kernel and a 
colored noise. We will examine the interacting scalar field theory that appears in the asymptotic 
limit of this non-Markovian process. To bring the reader to this new approach let us first consider 
the standard stochastic quantization of a free scalar field theory. 



3 Stochastic quantization for a free scalar theory 

Let us consider a neutral scalar field with a (Ay?^) self-interaction, defined in a d-dimensional 
Minkowski spacetime. The vacuum persistence functional is the generating functional of all vac- 
uum expectation vahic of time-ordered products of the theory. The Euclidean field theory can be 
obtained by analytic continuation to imaginary time supported by the positive energy condition 
for the relativistic field theory. In the Euclidean field theory, we have the Euclidean counterpart 
for the vacuum persistence functional, that is, the generating functional of complete Schwinger 
functions. Actually, the (A(^^)d Euchdean theory is defined by these Euchdean Green's functions. 
The Euchdean generating functional Z[h] is formally defined by the following functional integral: 

Z[h] = J [dip] exp (^-So + J d'^x h{x)^{x)^ , (7) 

where the action that usually describes a free scalar field is 



Sol 



d'x [l{dipr + lml^'ix)y (8) 



and the interacting part, defined by the non-Gaussian contribution, is 



Si[^]= I d'x^ p^x). (9) 
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In Eq.(7), [dip] is a translational invariant measure, formally given by [dif] = Ha; dip{x). The terms 
A and rriQ are respectively the bare coupling constant and the squared mass of the model. Finally, 
h{x) is a smooth function that we introduce to generate the Schwinger functions of the theory 
by functional derivatives. In the weak-coupling perturbative expansion, which is the conventional 
procedure, we perform a formal perturbative expansion with respect to the non-Gaussian terms of 
the action. As a consequence of this formal expansion, all the n-point unrenormalized Schwinger 
functions are expressed in a powers series of the bare coupling constant A [19] [20]. 

The aim of this section is to discuss the stochastic quantization of a free scalar field. It can 
be shown that it is equivalent to the usual path integral quantization. The starting point of the 
stochastic quantization to obtain the Euclidean field theory is a Markovian Langevin equation. 
Assume an Euclidean d-dimensional manifold, where we are choosing periodic boundary conditions 
for a scalar field and also a random noise. In other words, they are defined in a d-torus VL=T'^. 
To implement the stochastic quantization we supplement the scalar field ^p{x) and the random 
noise r/(;r) with an extra coordinate r, the Markov parameter, such that i~p{x) tp^r, x) and 
r]{x) r](T,x). Therefore, the fields and the random noise are defined in a domain: T'^ x 
Let us consider that this dynamical system is out of equilibrium, being described by the following 
equation of evolution: 

d S S 

— (^(t,x) = -j^Ux}=^{T,x) + r]{r,x), (10) 

where r is a Markov parameter, rjlr, x) is a random noise field and Sq is the usual free action 
defined in Eq.(8). For a free scalar field, the Langevin equation reads 

d 

— (/?(r, x) = -(-A + ml )ip{T, x) + r]{T, x), (11) 

where A is the d-dimensional Laplace operator. The Eq.(ll) describes a Ornstein-Uhlenbeck 
process and we are assuming the Einstein relations, that is: 

{rj{T,x)),^0, (12) 

and for the two-point correlation function associated with the random noise field 

{viT,x)viT',x')), =25ir-T')5%x-x'), (13) 

where ( ...)^ means stochastic averages. In a generic way, the stochastic average for any functional 
of (f given by F[(p] is defined by 



J[dr)]F[ip]exp 


-yd'^xjdrrj'ir, 


x) 


J [dr]] exp 




\jd'^xjdTrf{T, x) 





(14) 
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Let us define the retarded Green function for the diffusion problem that we call G{t — t',x — x'). 
The retarded Green function satisfies G{t — r', a; — x') = if r — r' < and also 



d_ 



G{r-T',x- x') = S'^{x- x')S{t - t') . 



(15) 



Using the retarded Green function and the initial condition (/^(r, x)\r=o — 0, the solution for Eq.(ll) 
reads 

ip(T, x)^ dr' / d'^x' G(t -t',x- x')ri(T', x'). (16) 
Jo Jn 

In the following we are interested in calculating the quantity {(p{T^x)(p{T',x'))r,. Using Eq.(12), 
Eq.(13) and Eq.(16), we have 

fmin{Ti,T2) r 

{ if in, Xi)if{T2, X2) )r,^2 dr' d^x'Gin - r', X^ - X') G{t2 - t', X2 - x'), (17) 

t/ " f2 

where min{Ti,T2) means the minimum of ri and T2. Using a Fourier representation, the two-point 
correlation function ( 9?(r, x)ip{r', x') ),i = D{t, x; t' , x') is given by 



D{r,x;r',x') = / d'^p^^ e-^"+<^^^-^'^ 



x—x') 



{p^ + ml) 

It is not difficult to show that Eq.(18) can be written as: 



(18) 



D 



1 



{t,x-t',x')- ^2 ^rfE ^(l 



.l)n 



niQ \ 2 



r 



l+n-l 



(19) 



where r —\ x — x' \ and is the modified Bessel function of order u. 

We can use the Fourier analysis to show that when the Markov parameters r and r' go to 
infinity we recover the standard Euclidean free field theory. Therefore let us define the Fourier 
transforms for the field and the noise given by ip{T, k) and ?7(r, k). We have respectively 



(27r)2 J 



and 



r)(T,k) = -^ [ d'xe-^'^^r^ir.x). 
(27ri2 J 



(20) 



(21) 
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Substituting Eq.(20) in Eq.(8), the free action for the scalar field in the {d+ l)-dimensional space 
writing in terms of the Fourier coefficients reads 

Soifik)] Uk)=v>{T,k) = ^J d'^k<f{T, k){k'^ + ml)ip{T, k). (22) 

Substituting Eq.(20) and Eq.(21) in Eq.(ll) we have that each Fourier coefficient satisfies a 
Langevin equation given by 

|-(^(r, k) = -{k' + mlMr, k) + v{r, k). (23) 

The solution for this equation reads 

(^(r, k) = exp (-{e + mDr) Lp(Q, k) + dr'exp (-(A;' + m^)(T - r')) r;(r', k). (24) 

Using the Einstein relation, we get that the Fourier coefficients for the random noise satisfies 

(?7(r,fc)), = (25) 

and 

( r]{T, k)rj{T', k') )r, = 2{2nY5{T - r') 5''{k + k'). (26) 

Before investigate the interacting field theory, let us calculate the Fourier representation for the 
two-point correlation function, i.e., {lp{t, k)Lp{T' , k'))^. Using Eq.(24), we obtain three contribu- 
tions to the scalar two-point correlation function. The first one is given by 

exp (-(A;' + ml)T + {k'^ + ml) r') (^(0, k)ip{Q, k'), (27) 

and decay to zero at long time. Let us assume that (^(r, A;)|,-=o = 0. There are also two crossed 
terms, each first order in the noise Fourier component given by 

i-t' 

2(/?(0,A;)exp(-(A;^ + m^)r) ds ex.p [-{k'^ + ml){r' - s)) r]{s,k'). (28) 

Since we are assuming the Einstein relations, i.e., {ri{T,x))jf = 0, on averaging on noise, these 
cross terms vanish. The final term is second-order in the noise Fourier component. Again, the 
solution subject to the initial condition (p{T, k)\r=o = can be used to give 

£ ds exp (-(fc^ + ml){T - s)) r){s, k) £ c/aexp (-(/c'^ + ml){r' - a)) r]{a, k'). (29) 



Again averaging on noise and using the Einstein relation given by Eq. (26) we have that this term 
becomes 

/• min{T,T') , . 

25\k + k!) rfsexp(-(A;2 + mo)(r + r'-2s)) . (30) 

Assuming that t — t' and using ( (^(r, k)if{T', k') ),,|r=T' = D{k, k'; r, t') we have 

D{k- r, r) = {2'Kf5\k + k')j^^^-^^ (l - exp (-2r{e + ml))) . (31) 

In the following, we are redefining the two-point correlation function as D{k] r, r) {2Tr)'^D{k; r, r). 
In the limit when r — > oo we recover the standard two-point function of the Euclidean free field 
theory. Before going to the next section, wc would like to mention the existence of more general 
Markovian Langevin equations. We can introduce a kernel defined in the d-torus. The kerneled 
Langevin equation reads: 

— (^(r,a;) = - j d%K{x,y)j^^\^^y)=^(^r,y) + v{r,x). (32) 

The second moment of the noise field will be modified to: 

( ri{T, x) riir', x') ), = 2(^(r - r') K{x, x'). (33) 
Choosing an appropriate kernel, it can be shown that all the above conclusions remain unchanged. 



4 Stochastic quantization for the {\ip^)d scalar theory 

In this section we will analyze the stochastic quantization for the (A<^^)d self-interaction scalar 
theory. In this case the Langevin equation reads 

d A 

— (/?(r, x) = -(-A + ml )^(r, x) - ^¥'^(^, x) + /^(t, x). (34) 

The two-point correlation function associated with the random field is given by the Einstein 
relation, while the other connected correlation functions vanish, i.e., 

{'n{Ti,Xx)r]{T2,X2)...r]{T2k-i-,X2k-i))n = 0, (35) 

and also 

{'n{Ti,Xi)...r]{T2k,X2k) )n = Y^ {v{n,Xi)r]{T2,X2) )r,{v{Tk,Xk)v{ri,xi) )^..., (36) 
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where the sum is to be taken over all the different ways in which the 2k labels can be divided 
into k parts, i.e., into k pairs. Performing Gaussian averages over the white random noise, it is 
possible to prove that 



lim ( (p{Ti, Xi)lf{T2, X2)...ip{Tn, X„) )^ = 



(37) 



where S{(^) = So{ip) + Si{(^) is the d-dimensional action. This result leads us to consider the 
Euclidean path integral measure a stationary distribution of a stochastic process. Note that the 
solution of the Langevin equation needs a given initial condition. As for example 



(38) 



Let us use the Langevin equation to perturbatively solve the interacting field theory. One way 
to handle the Eq.(34) is with the method of Green's functions. We defined the retarded Green 
function for the diffusion problem in the Eq.(15). Let us assume that the coupling constant is a 
small quantity. Therefore to solve the Langevin equation in the case of a interacting theory we 
use a perturbative series in A. Therefore we can write 



ifir, x) = (^(°) (r, x) + A(^(^) (r, x) + AV^'^ (t, x) + .. 



(39) 



Substituting the Eq.(39) in the Eq.(34), and if we equate terms of equal power in A, the resulting 
equations are 



d_ 



+ (-A^ + ml 



(P^^\t,x) = r)(T,x), 



+ (-A^ + ml ) 



d_ 



+ (-A^ + ml ) 



^(^nr,x) = -^(^(°)(r,x))\ 
^^'\r,x) ^ -^ {v^'\r,x))\('\T,x), 



and so on. Using the retarded Green function and assuming that 



|t=0 



(40) 

(41) 

(42) 
0, W q, the 



solution to the first equation given by Eq. (40) can be written formally as 

^(0) (r, x) = r dr' I G{t-t',x- x')r]{T', x') . (43) 
Jo Jn 

The second equation given by Eq.(41) can also be solved using the above result. We obtain 

iP^^\t,x) = -7^/ dri [ d'^xiG{T - Ti,x - xi) 
61 Jo Jn 



(44) 
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We have seen that we can generate all the tree diagrams with the noise field contributions. We 
can also consider the n-point correlation function {'^{Ti,xi)(^{T2,X2)---(p{Tn,Xn))rf Substituting 
the above results in the n-point correlation function, and taking the random averages over the white 
noise field using the Wick-decomposition property defined by Eq.(36) we generate the stochastic 
diagrams. Each of these stochastic diagrams has the form of a Feynman diagram, apart from the 
fact that we have to take into account that we are joining together two white random noise fields 
many times. 

As simple examples let us show how to derive the two-point function in the zeroth order 
( (p{ti, xi)(p{t2, X2) )!^\ and also the first order correction to the scalar two-point-function given by 
{(p{ri, Xi)(p{t2, X2) )^^^. Using the Eq.(16) and the Einstein relations we have 

rmin{Ti,T2) r 

{ ifiin, x,)ip{T2, X2) )(°) = 2 / dr' d^x' G{t^ - t\ - x') G{t2 - r', X2 - x'), (45) 

' JO Jn 

which is just Eq.(17). For the first order correction we get: 

{ip{X,MX2))^^^ = ~{J dXsJ dX,(G{X^ - X,)G{X2 - Xs) + 

G{X,-Xs)G{X2-X^)^7]{X^)(^l dX,G{X^-X,)r]{X,)^ )^. (46) 

where, for simplicity, we have introduced a compact notation: 

l^rlA^jdX, (47) 

and also y:>{T,x) = '^{X) and finally r]{T^x) = r]{X). 

The process can be repeated and therefore the stochastic quantization can be used as an alter- 
native approach to describe scalar quantum fields. We stress here that the stochastic quantization 
is based in the fact that although one starts with the system out of equilibrium, the Markovian 
Langevin equation forces it into equilibrium. Moreover, when the thermodynamic equilibrium 
is reached, the stochastic expectation values will coincide with the Schwinger functions of the 
Euchdean field theory. 
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Figure 1: Perturbative expansion for the scalar field where crosses denote noise fields. 




(a) (6) (c) 

Figure 2: The corrections up to one- loop to the two-point correlation function. 

We can represent Eq.(39) graphically as figure (1) (the random noise field is represented by a 
cross). Using this diagrammatical expansion, it is possible to show that the two-point correlation 
function up to one- loop level is given by figure (2), where we represent the retarded Green function 
by a line and the free two-point function by a crossed line. The rules to obtain the algebraic val- 
ues of the stochastic diagrams are similar to the usual Feynman rules. For instance the two-point 
function at one-loop level is given by 

{b)^-^S''iki + k2) J d''k£ dTG{ki;n - T)D{k;T,T)D{k2;T2,T). (48) 

(c) = -^6''{k, + k2) J d^kj^' dTG{k2;T2-r)D{k;T,T)D{k,;n,T). (49) 

A simple computation shows that we recover the correct equilibrium result at equal asymptotic 
Markov parameters (ri = r2 — >■ oo): 
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Obtaining the Schwinger functions in the asymptotic hmit does not guarantee that we gain a 
finite physical theory. The next step is to implement a suitable regularization scheme. A crucial 
point to find a satisfactory regularization scheme is to use one that preserves the symmetries of 
the original model. In the stochastic regularization method the symmetries of the physical theory 
is maintained. There are in general two different ways to implement the stochastic regularization. 
The first one is to start from a Langcvin equation with a memory kernel. It is known from the 
literature [21] that this method can at best only remove two degrees of divergence. Another 
possibihty is to smear only the noise field in the probability functional [22] [23] : 



{m)v 



JM^l-^Mexp 


-iJd'^xJdTj dr' ri{T, x)K]['ri{T', x) 


/ [dr]] exp 


-jJd'^xJdTj dr' r]{T, x)K^^rj{T', x) 





(51) 



where Kj^ is a memory kernel. In this case we change the Einstein relations of the noise field to: 

{r]{T,x)7j{T',x'))^ =2K^{t,t')5\x-x'). (52) 
The smearing function should be chosen such that, when A — > oo: 

lim K^{t-t') = 5{t-t'), (53) 

A— >oo 

recovering the usual theory. 

Since the Langevin equation is unaffected by the stochastic regularization, the physical field is 
the same as in the regularized case. However, the zeroth-order two-point correlation function is 
given by: 

D{k- r, r') = 

+ k') r ds r ds' G(k; r-s) G(k] t' - s')Ka{s - s') = 
Jo Jo 

2S'^{k + k') ds ds' exp[-{T + t' - s - s'){k^ + ml))KA{s - s'). (54) 



It is possible to prove that a necessary condition that the regularization function should 
satisfy in order to render the divergent loops finite is Kj^^t) 
kernels obeying this condition were proposed: 



K 



^ l^) = ^A^(A'Mrexp(-A^|r|) 



T=o= 0. The following series of 

(55) 
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For the case n = we obtain, for the free two-point correlation function: 



hm D{k;T,T) — ^ tt^^ ^- (56) 

.^oc ^ ' ' ^ (P + mg) {A^ + k^ + mD ^ ' 

Since the stochastic diagrams contains crossed hnes in its loops, we have that the ultraviolet 
divergences can be regularized choosing an appropriate n. Note that it is possible to use a different 
regulator of the type Kcj{t) = ^ar'^^^. This regulator scheme is quite similar to the analytic 
regularization of Bollini et al [24] and Speer [25]. The relation between the cutoff rcgularization 
and the analytic regularization procedure has been clarified by Kay [26] and also Svaiter and 
Svaiter [27] [28] [29] in a series of papers studying the Casimir effect. In the next section we will 
study a generalized Langevin equation with memory kernel and colored noise. 



5 Generalized Langevin equation with colored noise. 

It is well known that the Langevin equation given by Eq.(32) is only one particular choice in a 
large class of relaxation equations. To make sure that the generalized Langevin equation can also 
be used as a quantization tool, one must check that the process converges in the asymptotic limit 
and also that converges to the correct equilibrium distribution. Therefore, the aim of this section 
is to discuss the {Xf'^)d scalar field theory that appears if wc start from a Langevin equation with 
a memory kernel and also a colored noise. To proceed, let us introduce a Langevin equation with 
a memory kernel given by 

-^ip{T,x) = - ds Ka{t -s) ^ Ux)=^{s, x) + ri{T, x) , (57) 

where the stochastic random field 'r]{T, x) satisfies the modified Einstein relations. 

( r7(r, x) rjir', x') ), = 2Kj,{\t - t'\) 6''{x - x'). (58) 

In this case where Kj^{\t — t'\) has a width in the fictitious time, the description is Gaussian in 
spite of being non-Markovian. For the free scalar field we have that the generalized Langevin 
equation reads 

d r 

—Lp{T,x) = -J^ dsKA{T-s){-A + ml)(p{s,x) + 'r]{T,x). (59) 
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Using a Fourier representation for the scalar field and the random noise field we get 

^ip{r, k) = -(^2 + ml) £ ds Ka{t - s)ip{s, k) + r?(r, A;). (60) 
Following Ref. [16], we define the Laplace transform of the memory kernel: 

Kiz)^ r dTKAT)e-'\ (61) 

^0 

With the initial condition <^(r, A;)|t-=o = 0, the solution of the Eq.(60) reads: 

/•oo 

cpfr, k) = / dr' G(k,T-T') rjir', k) , (62) 

Jo 

where using the step function 9{t), the Green function G{k,T) is defined by: 

G{k,T) =Q{k,T)e{T). (63) 
The fl{k,T) function that appears in Eq.(63) is defined through its Laplace transform: 

Q{k, t) = (z + {k^ + ml)K{z)) \ (64) 



It is easy to see that, in the limit K\{t) — > S{t), we obtain the usual result for the Green function. 
From Eq.(62) and the modified Einstein relations, we get that the free scalar correlation function 
D{k; T, r') is given by: 

D{k;T,T') = 

^25\k + k') ds ds'G(k,T-s)G(k,T' -s')KJ\s-s' \) 
Jo Jo 

= 25'^{k + k') r ds r ds' n{k, r-s) n{k, r' - s') Kj^{\ s - s' |). (65) 
Jo Jo 

To proceed we have to write D{k; r, r') in a simplified way. Note that the double Laplace transform 
of the right hand side is given by: 

/* OO' /' OO /* 7" /* T ^ 

/ dre-'^ dre-'^l ds ds' Q{k,T - s) Q{k,T' - s') Kj^{\ s - s' \) ^ 
Jo Jo Jo Jo 

roo roo 

^n(k,z)n(k,z') ds ds'e-^' e-^' Ka(\ s- s' \). (66) 
JO Jo 
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Now, with simple manipulations, we get: 

rds'e-^'''e-"ms-s' |) = :^M±^. (67) 

Jo Jo z + z' 

Therefore, we get the identity: 

roc POO rT pt^ 

/ dre-'^ dr'e-'^ ds ds'n(k,T - s)n(k,T' - s') Ka(\ s - s' \) ^ 
Jo Jo Jo Jo 

= n(kMz'){!^^). (68) 

Remembering Eq.(64), we can show that: 

So, in parallel with result Eq.(67), we finally obtain a very simple expression for D{k;T,T') in 
terms of r). We have 

D{k; T, r') = 2 0^ I r - r' \) - n{k, r) Q{k, r')) . (70) 

Now, we need an expression for our memory kernel in order to investigate the convergence of 
Eq.(70). Prom Eq.(55), we will have, for n = 0: 

^A(r) = ^A'exp(-A2|r|). (71) 

Then, from Eq.(61), Eq.(64) and Eq.(71), and applying the inverse Laplace transform, we will 
obtain the following expression for the Q-function: 

"(*-)=(f-(f)-Kf))-(-4)- (-) 

where we have defined a real quantity (3 given by: 

(5 = A^2(F + mg)-A2. (73) 
Similarly, we will have, for the Green function: 
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From the results above, it is easy to see that the free two-point function will be given by: 



5'^{k + k') 



= 2 



(A;2 + ml) 



'A^ . / ,8(t-t') 



cos 



cxp 



A2 

2 



'A2 



sm 



/3r 
2 



cos 



2 



'A2 

.7 



sm 



cos 



r-r') 
exp(^-^(T + r')) 



(75) 



For T = r', we get: 



D{k;T,T)^2 



5\k + A;') 



(5t 



1 - I sin{ — ) + cos( — 



(A;2 + ml) 

So, in the limit r — > oo, we obtain the following result: 

5\k + A;') 



(5t 



exp(-A^T) . 



L»(A;;t,t) = 2 



(A;^ + ml) 



(76) 



(77) 



which, as we see, does not present an improved ultraviolet behavior. In fact, this equation is 
much similar to the usual equilibrium result, up to a constant. This is a first evidence that the 
ultraviolet divergences appearing in the perturbative series when we consider the self-interacting 
theory may not be regularized. Now, let us study the self-interaction (Av?'')^ scalar field theory 
within this non-Markovian approach. Now, the Langevin equation reads: 



d_ 



ip{T, x)^ - ds Ka{t - s) ((-A + ml)ip{s, x) + x)) + r]{T, x) 



(78) 



We can solve this equation by iteration as before. Then, after equating terms with equal powers 
in A, we get: 



^ / A 2 

— + (-A^ + ml 



<^(°'(t,x) = v{t,x), 



d_ 



+ (-A^ + ml ; 



V>^'\r,x) = £dsK4T-s) {v^"\s,x))\ 

^ ^) = - ^ r ^^^^ ~ ' ^^'^ ' 

and so on. The solutions of the first two equations can be written as: 

^(o)(r,x)= rdr' [ d'^x'G{r-T',x-x')ri(T',x'). 
Jo Jn 



— + {-A, + ml 



(79) 
(80) 
(81) 

(82) 
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and 

(p^^\r,x) = / dTi I d'^XxG{T - Ti,x - xi) [ dsK/^{Ti-s) 

o! JO Jn Jo 

Uo' In ^^^^ ~ ~ ^')^(^'' ^')) ' (S^) 

with the Green function given by Eq.(74). 

For the n-point correlation functions, the perturbation theory will be similar to the Markovian 
case, except that for each vertex in the stochastic diagram there is a memory kernel associated. 
For instance, the one-loop correction for the two-point function is given by {{b) and (c) given by 
figure (2)): 



(h)^--5''{h + k2) I d^'kj^' drj^ dsG(h,n-T)D(k;s,s)D{k2;T2,s)Kj,(\T-s\). (84) 

(c) = ~ + k2) I d'^k r dr r ds G{k2, T2 - T)D{k; s, s)D{ki; n, s)Ka{\ r - s \). (85) 
2 J Jo Jo 



From Eq.(74) and Eq.(75) we may split Eq.(84) in four parts as: 



where: 



(6) = - - S'^iki + k2) I d'^k £ dr £ ds [h + l2 + h + h) , (86) 

h= r dr r dsn{ki,Ti-T)n{k2,T2-s)n{k,0)KAi\T-s\), (87) 
Jo Jo 

h = r dr r ds n{ki, Ti - r) n{k2, T2) n{k2, s) n{k, 0)Ka{\t-s\), (88) 

JO JO 

Is= r dr r dsn{ki,Ti-T)n{k2,T2-s)n''{k,s)KA{\T-s\), (89) 

JO JO 

h= r dr r dsn{ki,Ti-T)n{k2,r2)n{k2,s)n\k,s)KA{\T-s\). (90) 

JO JO 

Using simple calculations, it is easy to show that the terms I2, I3 and J4 have a decaying oscillatory 
regime for Ti = r2 — > 00. So, we only have to consider the first term. It reads: 

/i = 

'A(ri-r) 



-exp(^--(r,+r2)jy^ dr(^-sm(^^ j + cos 



-(-f ) 1: - (f -(^) - -(^)) -(^^ ) • 
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After simple manipulations [30] and a tedious algebra, we obtain the final result, in the limit 
Ti = T2 ^ oo: 

where: 

f(\-n ... / (72A^ + 9/3?-/3|)A^ + /3|(/3?-/3|) U 2A^ \ 



and /3j = Ay2(k'f + mg) — A^; i = 1,2. So, we see that, although we obtain convergence in the 
asymptotic limit of the fictitious parameter r, we do not get a regularized theory. In fact, up to 
polynomials of A, the result is much similar to the usual equilibrium result, Eq.(50). Since the 
usual stochastic regularization requires the smearing of only the rj probability functional, leaving 
the Langevin equation alone, a natural question that arises is why a theory with a colored, internal 
noise does not lead to finite results in perturbation theory. Let us try to answer this question within 
a Fokker-Planck analysis. 



6 The Fokker-Planck approach 

Although in this paper we study the stochastic perturbation theory using the Langevin equa- 
tion approach, to understand our results it is more suitable to work within the Fokker-Planck 
formulation. As we know, correlation functions are introduced as averages over rj: 

Af J [dr]](^{Ti,Xi)Lp{T2,X2)---f{Tn,Xn)exp(^--Jd'^X J dr T]^ {t , x)j , (94) 

where J\f is given by 

jV-i = J [dri] exp(-^ J (fx j dTri\r, x)^ . (95) 

An alternative way to write this average is to introduce the probability density P[<f, r], which 
is defined as [8]: 

Pb,r] ^ J [dr^]expi^-y d'x J rfr' 77'(r', x)) U^i^iv) ' ^i^^v))- (96) 



where, for simplicity, we have absorbed the factor Af in the original functional measure. In terms 
of P[(fi,T], the equal time correlation functions will read: 



{(p{Ti,Xi)ip{T2,X2) ■ ■ ■(p{Tn,Xn))rj = / [d(p] (p{Ti, Xi)ip{T2, X2) ■ ■ ■ (p{Tn, Xn)P[^, t]. (97) 
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The probability density P[(p,T] satisfies the Fokker-Planck equation [31]: 
with the initial condition: 

P[(^,0]=n (99) 

y 

The stochastic quantization says that we shall have: 

where the limit is supposed to be taken "weakly" in the sense of the reference [8]. Again following 
the former reference, for the perturbation theory, we may split the action in two parts: 

S = So + XSi. (101) 

The Green functions of the theory are computed as power series in A. Similarly, we may expand 
the probability density as: 

oo 

P[^,T]^Y.Pk['P,r]. (102) 



So, the n-point correlation function in the kth order of perturbation theory will be: 

{(p{Ti,Xi)ip{T2,X2) ■ ■■(f{Tn,Xn) )rt = ^^ j [^^A ^{n, Xi)'^{t2, X2) ■ ■ ■ (f{Tn, Xn) Pk[p , t]. (103) 

The Fokker-Planck equation becomes: 
with the initial conditions: 

Po[</',0]=n %(?/))> (105) 

y 

Pfc[^,0] = 0, /c = l,2,--- (106) 
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r] = / W] [ .y r - r'] / ^ .'). (107) 



We may write a formal solution to the Fokker- Planck equation as: 

b 6Si 

where the Green functional of the free Fokker-Planck equation satisfies: 

with the boundary condition: 

Pob,^',T - r']|,=,; = n<^(<^(^) - <^'(^))- (109) 

X 

For our case, the above equations should be modified in order to take into account the presence 
of the memory kernel. For instance, the free Fokker-Planck equation will be: 

l-^PA^, r] =jAj A M(x -y,r) + j^^Ny ; r, x]) P„b, r], (110) 



M(x-y,r)=f dsKAr - s)GiT - s,x - y), (111) 
Jo 

N[lp ; r, = 5{x - y) I ds Ka{t - s)Sq \^(x)=^{s,x) ■ (112) 



where: 
and 



Now, we can easily see why we required earlier that defined in Eq.(73), should be real. Re- 
membering the familiar transformation, P[(/9,r] = ■?/'[(/?, r] exp(^—^^^, with S = M'-^N (M — 
M{x — y,T), N = N[^p ; r, x]), we can recast Eq.(llO) into the Schroedinger type equation: 

^ = -2ni^, (113) 

where Ti is the Fokker-Planck Hamiltonian: 

n=]^j d^xj d''yQ{x)M{x-y,T)Q{y), (114) 

where: 

Q{x) = Il{x) + -[fl{x),s], (115) 
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and n(a;) = —ij;^- We stress here that this derivation is purely formaL Now, let us study the 
function M{x — y,T). In the Fourier space, it is given by: 



4 exp( — - 



2A^cos(:-) + (:^-^Jsin^2 



(116) 



where we have used Eq.(71) and Eq.(74). So, we see that, in the limit r — > oo, the M-function will 
be positive definite which, according to Eq.(114), guarantees that the Fokker-Planck Hamiltonian 
be positive definite as well. In other words, the real part of its eigenvalues will be positive. But 
what happens if /3 is a purely imaginary number? As for the Green function, we will have: 

G{k,r) = f^sinhf^) +coshf^')^ expf-r4^V(r). (117) 



P \ 2 J \ 2 J J "^K 2 

So, from Eq.(116) and Eq.(117), it is easy to see that M{k,T) may not be positive definite and, 

therefore, neither will be Ti. Then, there is no guarantee that the real part of the eigenvalues of 
the latter will be positive. In other words, we cannot assure that, in the limit r — >■ oo, the system 
will reach its ground state, i.e., that it converges to an equilibrium. 
Now, let us study the solution to Eq.(llO). It will be of the form: 

Pokr] = f[dv]l[Kip{y)-<p{r,y)) 

J y 

d'^x J dr" J dT'ri{r",x)Kx\r" -T')r]{T',x)^, (118) 
where (p{T,y) satisfies Eq.(62). We can rewrite this last equation in the form: 
Po[v,r] = / [dri] I [dC]e^v(i J d^x C{^{x) - ^{t,x))^ 

exp(^-^ J d'^x J dr" J dr' r}{r", x) Ki\t" - r') r]{T', x)^ . (119) 

Using Eq.(62) and doing the r], ^ functional integrations, we obtain: 

Po[ip,T]= No^exp(^-^J d'^x J d'^x'<^{x) D-^ {T,x;r,x')<^{x')^, (120) 
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where the D-function is just the free two-point correlation function and Nq^ is a normahzation 
factor. In momentum space, we get, for r — oo: 

Po[</^,r]=iVo-iexp(-^y d'k ipik){e + mlM-k)). (121) 

As we can see, this result resembles the usual equilibrium result, up to a constant, in contrast with 
the result obtained by Breit et al [22]. Now, proceeding with similar steps (see reference [8]), it is 
possible to show that the Green functional for the free Fokker-Planck equation is given by (Ag ^ 
is a normalization factor): 

Vo[<p,ip',T] = Ao-^exp(^-i J d''kipi{k)Dij{r)ipj{-k)'^, (122) 

where: 

ip. = i^ik) if'ik)), (123) 
and the elements of the matrix Dij are given by: 

D(r)^( D(k;r,r) -D(k;r,r)G(k,r) \ 

' \-D(k;r,r)G(k,T) D{k;T,T) G(k,r) G(k,T) J ' ^ ' 

It is easy to verify that, in the limit r — > oo, Fq and Vq will satisfy, up to constants, similar 
relations as obtained by Floratos and Iliopoulos [8] . This behavior can be understood if we notice 
the similar structure between the Green function of the usual Parisi-Wu scheme and our non- 
Markovian approach. The same is true for both free two-point correlation functions. However, 
for massless scalar theories, those estimations in [8] decay as inverse power of r. In our approach, 
even in this massless situation, an exponential behavior is found. Therefore, in the limit r — > oo, 
we get an improved convergence. 

Prom all of these results we expect to obtain probability functions in the perturbation theory 
that are very close to those usual equilibrium probability densities. Therefore the theory will not 
be regularized. As remarked on [16], nothing is lost by using a non-Markovian description in place 
of a Markov description, as long as it is Gaussian. 
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7 Conclusions 



It is well known that the Langevin equation given by Eq.(32) is only one particular choice in 

a large class of relaxation equations. The aim of this paper was to investigate if the Parisi-Wu 
quantization method can be extended assuming a Langevin equation with a memory kernel with 
the modified Einstein relations. Therefore, in this paper we discussed the stochastic quantization 
for massive self-interacting scalar field, first assuming a Markovian Langevin equation. Then, we 
studied the (A(^^) theory at the one-loop level introducing a non-Markovian Langevin equation and 
examine the interacting field theory that appears in the asymptotic limit of the non-Markovian 
process. 

To make sure that the first modification can be used, one must first check that the system 
evolves to the equilibrium in the asymptotic limit. Second we have to show that converges to 
the correct equilibrium distribution. We proved that although the system evolves to equilibrium, 
we obtain a non-regularized theory. In contrast with dimensional regularization [32] [33] [34] [35] 
and also other regularization procedures, a remarkable property of the stochastic regularization is 
that it preserves all the symmetries of the original theory. With our results, it is easy to see that 
the system described by this non-Markovian Langevin equation converges, but the virtues of the 
stochastic regularization that we hoped to appear in this framework are lost. 

A natural continuation of this paper is to discuss the stochastic quantization of bosonic and 
femionic fields in general Riemannian spaces. A different application of the stochastic quantization 
is to discuss interacting field theory in the presence of macroscopic structures. It is well known 
that performing the weak-coupling perturbative expansion, to renormalize the interacting field 
theory, we have to introduce not only the usual bulk counterterms, but also surface counterterms. 
This can be done at the one-loop level at zero and finite temperature [36] [37] [38]. To extend 
the calculation to high-order loops will appear overlapping divergences. A natural question is how 
to implement the stochastic quantization in systems with macroscopic structures. This subject is 
under investigation by the authors. 
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